7r — Transition in Superconductor-Ferromagnetic-Superconductor Junctions 
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Superconductor-Ferromagnetic-Superconductor 
(SFS) Josephson junctions are known to exhibit a 
transition between n and states. In this note we find 
the 7T — phase diagram of an SFS junction depend- 
ing on the transparency of an intermediate insulat- 
ing layer (I). We show that in general, the Josephson 
critical current is nonzero at the n — transition tem- 
perature. Contributions to the current from the two 
spin channels nearly compensate each other and the 
first harmonic of the Josephson current as a function 
of phase difference is suppressed. However, higher 
harmonics give a nonzero contribution to the super- 
current. 

PACS numbers: 74.50.+r, 74.80.-g, 75.70.-i 

In the last years, many interesting phenomena were 
investigated in Superconductor (S) - Ferromagnetic (F) 
- Superconductor Josephson contacts. One of the most 
striking effects is the so called 7r-state of SFS junctions 
[Q-Q in which the equilibrium ground state is character- 
ized by an intrinsic phase difference of ir between the two 
superconductors. Investigations of n junctions have not 
only academic interest, e.g., in a solid-state imple- 
mentation of a quantum bit was proposed based on a 
superconducting loop with and n Josephson junctions. 

The existence of the 7r-state in an SFS junction was 
recently experimentally demonstrated by the group of 
Ryazanov In this experiment, the temperature de- 
pendence of the critical current was measured. At a cer- 
tain temperature the critical current was found to drop 
almost to zero, this has been interpreted as the transition 
from the to the 7r state. The transition temperature T^q 
was shown to exhibit a strong dependence on the concen- 
tration of ferromagnetic impurities, i.e., on the exchange 
field E cx in the ferromagnetic film. 

In this note, we present a theory of the ir — transition 
in short SFS junctions. Our goal is to understand which 
parameters (exchange field, temperature,...) stabilize a 
7r state and how the phase diagram looks like. We inves- 
tigate the current-phase relation and the critical current 
near the transition to the 7r-state. Most importantly, we 
find that in general, the critical current is not zero at 
T^o, and it may not even reach a local minimum. The 
identification of the critical current drop and the tt — 
transition is only possible if the current is given by the 
standard Josephson expression I(ip) oc sin(y), which is 



valid for the limiting case of tunnel barriers only. Even 
if the main contribution to the current is of this form, 
the higher harmonics contribution I(ip) oc sin(2c/?) would 
not vanish at T^o- Consequently, I c ^ at the transition 
point. 

We consider SFS junctions in the "short" limit de- 
fined by H/t ^> A(T = 0); here, r is the characteris- 
tic time needed for an electron to propagate between the 
superconductors. In this case, we can employ a pow- 
erful scattering formalism which allows one to ex- 
press the energies of Andreev states in the junction in 
terms of the transmission amplitudes of the junction in 
the normal state. These Andreev states give the main 
contribution to the phase-dependent energy of the junc- 
tion, therefore I(ip) can be calculated. Any junction is 
characterized by a set of "transport channels" labelled 
by n = 0,1,... ,N, each channel is characterized by 
the transmission coefficient D n . If one disregards fer- 
romagnetism, the Andreev levels are degenerate with re- 
spect to the spin index a. Their energies are given by 
E na = ±A(l-A l sin 2 (<^/2)) 1 /2. 
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FIG. 1. Junction configuration. The current flows from 
one superconductor (S) to the other through the ferromag- 
netic (F) layer (width d) which a scattering region denoted 
by I. The exchange field is supposed to be parallel to the SF 
boundaries. 

We generalize the scattering approach to cover SFS 
junctions. In this case, the phase of the transmission 
amplitudes also becomes important. To see this, we in- 
troduce the parameter 7: cos(-f(ip)) = 1 — 2D n sm 2 (ip/2). 



1 



If we assume that ferromagnetism does not change the 
transport channels, the energies of the Andreev states 
become 



(1) 



where $> n ,a is the phase of the transmission amplitude 
for an electron with spin a in the channel n. Thus we 
observe that the different phase shifts for different spin 
directions result in a spin-dependent energy shift of the 
Andreev states. To specify the model further, we con- 
sider the layout shown in Fig. It consists of two bulk 
superconductors, a ferromagnetic layer with exchange en- 
ergy E cx <C Ep, and a scattering region denoted by 
I. We assume that the F-layer is ballistic and that the 
order parameter A is constant in the superconductors: 
A(x) = Ae ±v / 2 , and A (at) = in F. 

We believe that the model considered is quite a gen- 
eral one. It is applicable to quasiballistic SFS mul- 
tilayers (recently a quasiballistic SF junction was pre- 
pared by Kontos et al. ||) with either specular or dis- 
ordered interfaces ||; also to Josephson junctions where 
electrons tunnel through small ferromagnetic nanoparti- 
cles |l(],[ll]]. We shall restrict ourself to the case, when 
the width of the scattering region is much smaller than 
the width d of the junction. The transport channels 
can be associated with different incident angles. Then 
3> W)CT — &n.-a — cnr(2E CK d/TrhvF)ln = CT7r0Z„, where 
l n > 1 is the length of a quasiparticle path between the 
superconductors divided by d, see Fig. |. Formula (0) 
reproduces the energy spectrum obtained in the limiting 
case D = 1 in 

The contribution to the free energy of the junction 
which depends on ip is given by 



cosh 



2T 



(2) 



The continuous spectrum is neglected in (En; one can 
easily check that it gives a (^-independent contribution 
to the free energy. The summation over the channels 
n can be evaluated by converting the sum to an inte- 
gral: = J dlp(l), where p(l) = Y, n $( l ~ l n) and 
J p(l)dl — N, the number of channels. If there is only 
one channel in the junction, the weight function p de- 
fined above reduces to 8(1 — 1). If, on the other hand, 
the number of channels N is much bigger than unity, 
p(l) = 2N/l 3 9(l - 1). (We assumed D to be indepen- 
dent on n.) A similar distribution of / can be found for 
SFS junctions with disordered boundaries, see 0. At 
some points of these notes, we shall use the distribution 
p(l) = N8(l — 1), since it allows us proceed analytically, 
and the results obtained with it are qualitatively the same 
as with the other distributions. We shall refer to this dis- 
tribution as the ^-distribution. (When p(l) — N8(l — 1), 



our parameter is closely related to spin-mixing angle 
introduced in PH| .) 

Which exchange field in F is sufficient to ensure that 
the SFS junction can be put into a 7r-phase by changing 
the temperature? The 7r-state is the result of the ferro- 
magnetic exchange field in the F-layer. If it is too small, 
then the junction will remain in the 0-phase at all tem- 
peratures. We show below which values of exchange field 
and temperature guarantee that the junction will be in 
the 7r-phasc. 

In an equilibrium situation with zero current, the tem- 
perature Ttto separating the 7r and phases is determined 
from the condition that the free energy f2 reaches its min- 
imum at (p = 2im and at ip = T7 + 2irn, n = 0, ±1, . . . (the 
free energy of an ordinary junction has a global minimum 
at ip — 2ixn). The numerical solution of this equation 
for T v o(D) is shown in Fig. ||. Here and below, we use 
the approximation A(T)/A(0) = tanh(1.74 v /T c /T - 1) 
in doing numerical calculations. 
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FIG. 2. The temperature of the it — transition at zero cur- 
rent versus the dimensionless exchange field O = 2E cyi d/iihvF 
at different values of junction transparencies D. Only trajec- 
tories with / = 1 are taken into account which is justified 
if either one channel or many channels are present. Inset: 
phase diagram of the junction at D — 0.1. The gray regions 
correspond to the 7r-phase, the white regions to the 0-phase. 

If D = 1, the 7r-phase can exist only in the domain 
2n + 1/2 < 6 < 3/2 + 2n, n = Q,±l,.... At finite D, 
there are regions of in which either the 7r-phase or the 
0-phase is stable for all temperatures, I C (T) has no cusps 
in these regions. For — > 1/2+n, T n0 — > T c for arbitrary 
transparency. 

There are regions in the phase diagram where f2 has 
two minima, it — 2im and ip — n(2n + 1), n = 0, ±1, .... 
We shall consider these regions below. 

The (T, 0) phase diagram of the junction is depicted 
in the inset of the Fig. |2[ The diagram is periodic in 
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with period 2tt. It follows from the graph that large value 
of the exchange field = 2E C xd/irhvF do not guarantee 
that the SFS junction is a 7r-junction. 

Evidence of the existence of a tt phase in SFS junc- 
tions was experimentally demonstrated by the group of 
Ryazanov pj. The experimental curves I C (T) showed 
cusps at a certain temperature (which we shall denote by 
T^o); the critical current at the cusp was close to zero. 
There are qualitative arguments in p| that the cusp cor- 
responds to transition of the junction to the tt state and 
I c = at the cusp. We agree with the first statement, but 
disagree with the second. In our opinion, there is no qual- 
itative argument for the critical current to be zero at the 
temperature of the cusp. Our model gives qualitatively 
similar curves I C (T) to those presented by Ryazanov et 
al, but I c 7- at the cusp, where the junction undergoes 
transition between and n states. This will be discussed 
below in more detail. 

The Josephson current I carried by the Andreev levels 
(|l|) can be found from the free energy (||) using the re- 
lation I(tp) = ^d v fl((p). Using d v "f = D sin ((f)/ sin (7), 
we obtain 



w \ v — v — ^ 2e D sin((p) . 



h sin(7) 



7 + O"7T0Z„ 



(3) 



tanh 



— cos 

2T \ 2 



If = and D = 1, Eq. (Q) reduces to the usual for- 
mula for the Josephson current in a short ballistic SNS 
junction, leading to the critical current I c = NeA/h at 

If D <§C 1, we can proceed analytically in the calcula- 
tion of I C (T). Then 7 ps 2a/D| sin((p/2)| < 1 and we can 
expand the Josephson current (||) in 7. This leads to 



eAD 
2h 



(4) 



where 



/(e,T) 



2T 



„ m f / ^, / n\ , /cos(tt0Z/2) A 
dlp{l) < cos (nQl/2) tanh — v 

Asin 2 (vreZ/2) 
2Tcosh 2 (cos (ttOZ/2) A/2T) 



(5) 



If g > in (|^), the junction is in the ordinary state. In 
the opposite case, the junction is in the tt state. 

Solving g(Q,T) = gives us the transition temper- 
ature TVo- For p(l) oc 5(1 - 1), g(Q,T) = can be 
solved only in the domain 2n + 3/2 < O < 1/2 + 2n, 
n = 0, ±1, .... As a consequence, the 7r-phase exists 
only in these domains. If — > 1/2 + n then T — > T c ; 
0^1 + 2n leads to T -> 0. 

In the region |T — T^ol ~ 2? T c , the current is no longer 
given by Eq. (|J). Higher harmonics in tp have to be taken 
into account. Since the n-th harmonic is proportional to 



D n sm(nip) (follows from (^)) and D <C 1, the second 
harmonic gives the main contribution to the Josephson 
current: I((p) oc D 2 sm(2ip). That means the critical 
current is not zero at the temperature of n — transition, 
I c oc D 2 . Near T^o the currents of the spin channels a = 
±1 flow in opposite directions and nearly compensate 
each other; therefore I c is suppressed. 

Near the critical temperature T c , we also can proceed 
analytically. Then we obtain from (j^) for the Josephson 
current: 



2Th 



sm(ip) / dlp{l)Dsa\(TTQl) 



(6) 



Using the ^-distribution of the trajectory lengths, we find 
that the junction is in the 7r-state when 1 + 2n < 9 < 
2 + 2n, n = 0,±l, .... 

Figure ^ shows the typical dependence of the critical 
current on temperature. 
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FIG. 3. Dependence of the critical current on temperature 
for <d = 0.7 at different transparencies D. Cusps in the I C (T) 
curves indicate the 7r — transition. I c is always nonzero at 
the cusp, see the inset. The plot of the critical current for 
D — 0.22 has two cusps. At low temperature and I < I c , 
the junction with D = 0.22 is in the 7r-state, at intermediate 
temperatures (between two cusps) - in the 0-state, and for 
temperatures near T c - again in 7r-state. 

The critical current in the figure is normalized to the 
critical current I c q — N(eA/h)(l — yl — D) at zero tem- 
perature and zero exchange field. The plot corresponding 
to D = 0.22 exhibits two cusps. When the temperature 
is low the junction is in the 7r-state; for intermediate tem- 
peratures (between the two cusps) the junction will be in 
the 0-state, for T <T C - again in the 7r-state. (There is 
a schematic plot in ffl] where I C (T) of a SFS junction has 
many cusps with nonzero critical current in the cusps. 
However, Ref. |jj does not provide an explanation of this 
fact). The critical current for D = 0.3, 0.4, 0.7 in Fig is 
greater than the critical current I c q corresponding to zero 
temperature and zero 0. The exchange field enhances the 
Josephson current in the SFS junction |14|. 
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Below we shall investigate the ir — transition when 
a d.c. current I < I c is injected into the junction. We 
shall pay attention to the regime when / is smaller than 
I c at the cusp. Suppose that the temperature is changed 
at fixed i\ Then at the 7r — transition temperature, the 
phase across the junction will jump approximately by 7r. 

There are several solutions f(T) of the equation / = 
I{tp), where I{(p) is given by (|). We will assume that 
the damping is large, such that the phase is stabilized 
at one of the minima of the Gibbs energy Zz(I,T,tp) — 
f2(T, (p) — (plh/2e [jl5|. The phase values corresponding 
to the minima of the Gibbs energy are depicted by solid 
lines in Fig. |]b. If the temperature is increased from 
T = 0, the phase will continuously change with T until 
T reaches the dark gray region in Fig^ja where 5 has two 
minima of tp in [— tt, tt]. Here the phase will choose one of 
the minima depending on the dynamics of the junction, 
which depends on the properties of the external circuit. 
Outside of this region at higher temperatures, the phase 
will also follow continuously adiabatic changes of T. 



In conclusion, we investigated the phase transition be- 
tween the 7r and phases in a ballistic SFS junction with 
a scatterer in the F layer. We calculated the (T, 8) and 
(J, T) phase diagrams of the junction. It was shown that 
there is no reason for the critical current to be zero at the 
transition temperature T n Q. The currents of the two spin 
channels nearly compensate each other at T^o, and the 
current scales as D 2 sin(2</j), D <C 1 instead of Dsin(ip), 
as it does far from the transition. 
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and M. V. Feigel'man for stimulating discussions and use- 
ful comments on the manuscript. Our research was sup- 
ported by the Swiss National Foundation; N. M. C. is 
also supported by the RFBR, project No. 01-02-06230, 
by Forschungszentrum Julich (Landau Scholarship), and 
by the Netherlands Organization for Scientific Research 
(NWO). 
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FIG. 4. (a) The phase diagram of the junction for D = 0.1, 
O = 0.7. The light gray region corresponds to the 7r-phase, 
the white to the 0-phase. The Gibbs potential has two min- 
ima of tp in [— tv, 7r] in the dark gray region. The critical 
current (solid thick line) is the upper boundary of the phase 
diagram, (b) Temperature dependence of the phases corre- 
sponding to the d.c. current / = 0.07/ c o (dashed line parallel 
to the temperature axes in (a)). The thick solid line repre- 
sents the stable solution of the equation I = I(ip) (minimum 
of the Gibbs energy), the dashed curve exhibits the unstable 
solution (maximum of the Gibbs energy). The equilibrium 
transition temperature (I = 0) corresponds to T/T c = 0.2. 
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